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Urban heating in northern China accounts for 40% of total building energy usage. In central heating systems, heat 
is often transferred from heat source to users by the heat network where several heat exchangers are installed at 
heat source, substations and terminals respectively. For given overall heating capacity and heat source tempera- 
ture, increasing the terminal fluid temperature is an effective way to improve the thermal performance of such 
cascade heat exchange network for energy saving. In this paper, the mathematical optimization model of the cas- 
cade heat exchange network with three-stage heat exchangers in series is established. Aim at maximizing the cold 
fluid temperature for given hot fluid temperature and overall heating capacity, the optimal heat exchange area dis- 
tribution and the medium fluids’ flow rates are determined through inverse problem and variation method. The 
preliminary results show that the heat exchange areas should be distributed equally for each heat exchanger. It al- 
so indicates that in order to improve the thermal performance of the whole system, more heat exchange areas 


should be allocated to the heat exchanger where flow rate difference between two fluids is relatively small. This 


work is important for guiding the optimization design of practical cascade heating systems. 
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Introduction 


Heat exchangers and corresponding heat exchange 
networks are playing an important role in various engi- 
neering fields involving energy production, transfer, 
conversion and terminal utilization!’ *!. With the rapid 
industrialization and modernization over the recent two 
decades, the global total energy consumption has grown 
by 49 %"!. Therein, buildings account for more than 30% 
of the total energy consumption and the percentage keeps 
increasing due to the rapid urbanization"! During the 
same period, building energy consumption in China has 
increased with an average annual rate of over 10%!°!. 
According to the latest statistics analysis, urban heating 
in northern China accounts for 40% of total building 
energy usage'®!, Hence, improving the energy efficiency 
of heating systems is of high significance in building 
energy saving. Due to high energy and economic effi- 
ciency, low pollution emissions and high energy supply 
safety and reliability, central heating (e.g., combined 


heating and power (CHP) system) constitutes the main 
system form for space heating in northern China!””!, 
Cascade heat exchange networks containing several heat 
exchangers in series are widely used in central heating 
systems''”!, Generally speaking, heat is often transferred 
from heat source to users by two water circulating net- 
works (i.e., primary and secondary network), and three 
kinds of heat exchangers installed at heat source, substa- 
tions and terminal users respectively” ”. 

Many researchers dedicated to improve the thermal 
performance of central heating systems for energy saving. 
Hasan et al. used low temperature water for space 
heating via utilizing the floor radiator as the terminal 
devices. Zsebik et al.''*! found that a reasonable arrange- 
ment and configuration of heat exchangers in substations 
could effectively reduce the heat transfer loss. Khan et al. 
l4 and Barbieri et al." introduced thermal energy sto- 
rage into CHP central heating systems to bridge the 
timely gap between heat supply and demand. Fu et al.!’*! 
proposed the concept of absorption heat exchanger by 


integrating absorption heat pump with traditional coun- 
ter-flow heat exchanger, and applied the new system to 
practical central heating systems to take place of tradi- 
tional heat exchangers at heat source and substations re- 
spectively. The demonstrative project in northern cities in 
China showed that the new system could dramatically 
increase the heating capacity of the whole system for the 
same heat source amount. Based on absorption heat ex- 
changer, Wang et al.''”! put forward the absorption tem- 
perature transformer that could make a large temperature 
lift/drop through multi-stage heat exchange processes to 
further enhance the heating system performance. Then 
Zhang et al.''*! built the ideal central heating system con- 
taining heat engine, heat pump and heat exchanger, and 
presented that absorption heat exchanger served as a typ- 
ical form of the so-called heat adaptor. For the heat 
transfer enhancement and thermal performance im- 
provement, Shah et al.''®! established the relationship 
between entropy generation and heat exchange effec- 
tiveness for complex flow arrangement to increase heat 
transfer. However, Guo et al.”°! stated that the entro- 
py-based analysis method is only feasible for heat-work 
conversion processes and devices, but not applicable to 
heat transfer processes and corresponding equipment. So 
they proposed a new concept, entransy, for heat transfer 
processes analysis and optimization’, Liu et al.!7! 
compared entropy analysis and entransy analysis in heat 
exchanger optimization and found that it is entransy dis- 
sipation variation that always kept consistence with the 
optimal heat exchanger performance, instead of entropy 
generation. Then entransy-based analysis method has 
already been used by more researchers for heat exchang- 
er and heat transfer involved system optimization”*”*!, 
Nonetheless, existing researches almost focused on the 
optimization of specific heat transfer process or single 
heat exchange unit. Central heating system is often a 
comprehensive heat exchange network containing mul- 
ti-stage heat exchangers in series, so that not only the 
single heat transfer unit but also the system arrangement 
and configuration have a great impact on the system 
energy efficiency. Thus the optimization design of such 
cascade heat exchange network is of high significance in 
improving the thermal performance of the whole system. 
In typical central heating systems, heat is often trans- 
ferred from heat source to users by three-stage heat ex- 
changers and two medium fluids circulation networks. So 
how to determine the optimal heat exchange area distri- 
bution and medium fluid flows is an important problem. 
In this paper, aim at maximizing the cold fluid tempera- 
ture for given overall heating capacity and hot fluid tem- 
peratures, the simplified optimization model of the cas- 
cade heat exchange network is established and the op- 
timal heat exchange areas and the medium fluids’ flow 
rates are obtained through inverse problem and variation 


method. This work is important for guiding the optimiza- 
tion design of practical cascade heat exchange systems. 


Method 


Problem Description 


The typical process of a central heat system is shown 
in Fig. 1. Different kinds of heat sources, such as coal/gas 
boiler, extracted steam from thermal power plant, low 
grade industrial waste heat and so on, are utilized to pro- 
duce high temperature water supply by heat exchangers 
near the heat source. Then the hot water is delivered to 
the heating zone through the primary heat network and 
the heat is transferred to the water in the secondary net- 
work with temperature decrease by the heat exchangers 
at district substations. And the relatively low temperature 
water is pumped to the users and the heat is transferred 
again to the indoor air for space heating by the terminal 
heat exchangers such as radiators and fan coils. 
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Fig. 1 Schematic diagram of central heating system. 


So a typical central heating system is indeed a cascade 
heat exchange network with three-stage heat exchangers 
in series (i.e., heat source, substation and users) and two 
medium fluid circulations (i.e., primary and secondary 
heat network). For analysis simplification, as Fig. 2 
shows, such a central heating system can be regarded as a 
cascade heat exchange network, containing three heat 
exchangers and two medium fluids. Firstly, high temper- 
ature heat source (Thi, Thos Cy) is utilized to produce hot 
water in the primary network (7), T>e C2). Secondly, 
heat is transferred from the primary network to the sec- 
ondary one (Tim Tic Cı) at the substations. Finally, the 
water flows into the terminal devices to heat the indoor 
air (Toi, Te,os Ce). 
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Fig. 2 Cascade heat exchange network with three heat ex- 
changers and two medium fluids. 


For space heating, the main objective is to increase the 
terminal indoor air temperature (Toi, Te,o) for given over- 
all heating capacity (Q) and heat source temperature (Thi 
Tho). Therefore, the thermal performance of the whole 
system highly depends on the cascade heat exchange 
network design. When the overall heat exchange areas 
are given ((UA),+(UA),,+(UA),=(UA)), how to allocate 
them for the three heat exchangers ((UA),, (UA)m, (UA)¢) 
If take economical factor into consideration, for instance, 
the middle heat exchanger is more expensive than the 
others so that its heat exchange area needs to keep a cer- 
tain value ((UA)m<(UA),, (UA).), what is the difference 
in heat exchange areas allocation and how to determine 
the medium fluids flow rates? To answer these questions, 
the optimization problem of the optimal heat exchange 
areas distribution and medium fluids flow rates determi- 
nation for such a system needs to be solved through the 
following approach. 


Inverse Problem and Variation Method 


Compared to forward problem approach, inverse 
problem method shows obvious advances in solving 
some optimization problems in different engineering 
fields'’® 7”!, Based upon that, Zhang et al."! proposed a 
new approach, by combining the inverse problem method 
and Lagrange variation method. This approach shows 
special advantages in solving building energy and envi- 
ronment optimization problems. The essence of the in- 
verse problem and variation method is the “least action 
principle”, and the specific procedure can be roughly 
divided into two steps: (1) construct the Lagrange func- 
tion (F(x;)), which is also regarded as the “action” of the 
given problem, according to the optimization objective 
(f(x;)) and constraint conditions (g,(x;)); (2) make the “ac- 
tion” equal its extreme value through variation method 
(5F(x)=0), to obtain the undetermined parameters (x;). 


The mathematical expression of the inverse problem and 
variation method are as follows: 
Undetermined parameters: x,,i=1---n; 


Objective: min/ max f(x); 
Constraint condition: g;(x;)=0,j=1---m. 
In order to converse the constraint optimization problem 


to an unconstraint one, the Lagrange function can be built: 


F(x) = fat’ a) (1) 
J 


where A; is the undetermined Lagrange coefficient. Thus 
the equivalent objective of the optimization problem can 
be obtained: 
min/ max F'(x;) <> min/ max f(x;) (2) 
Then regard the constructed F(x;) as the “action” and 
make it equal to the extreme value through variation me- 
thod: 
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Finally, by solving the above equations, the desirable 
parameters can be obtained under the given objective and 
constraint conditions. 


Optimization Model and Solution 


For the aforementioned cascade heat exchange net- 
work optimization problem (Fig. 2), the overall heating 
capacity (Q), heat source temperatures (This Tho), flow 
rates of both hot and cold fluids (Ch, Ce) are all known 
parameters. The objective is to determine the optimal 
heat exchange area distribution of the three heat ex- 
changers ((UA),, (UA), (UA).) and the flow rates of the 
two medium fluids (C4, C2), in order to obtain the highest 
cold fluid temperatures (Toi, Teo). According to the in- 
verse problem and variation method, the optimization 
model of such a problem can be established as follows: 

Known: QO, This Thos Ch Ce 

Unknown: (UA), (UA) im, (UA)e, Ci, Cr} 

Objective: max T,;; 

Constraint condition: (UA),+(UA)_+(UA)=UA. 

Based on the optimization objective and constraint 
condition, the Lagrange formula for this problem can be 
expressed by: 

F =T,,+A{(UA), +(UA) +A), -UA] (A) 

According to energy conservation for each heat ex- 
changer, there is 
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From Eqs. (4) and (5), by making 5F=0, it can be ob- 
tained that 
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Then the optimal heat exchange areas and medium 
fluid flow rates can be calculated out: 
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The aforementioned problem solution is under the ideal 
situation that there is no more restrictions, such as eco- 
nomic constraint conditions, on those heat exchangers. If 
take economic factor into account, for instance, suppose 
that the middle heat exchanger at the substation is much 
more expensive than the others, so that (UA), should keep 
a small and certain value ((UA)m<(UA)p, (UA),). Then the 
updated optimization model can be expressed by 

Known: Q, This Thos Chs Cos (UA) mn’ 

Unknown: (UA), (UA),’, Cy’, C2’; 

Objective: max Toi’; 

Constraint condition: (UA),’+(UA),’=UA. 

Similar to the previous case, the Lagrange formula can 
be established: 


(7) 
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F =T, ,'+ (UA), '+ (UA).'-(UA)'] (8) 
Through the same method, it can be obtained that 
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So after considering the extra restriction, the optimal 
parameters vary, compared to the previous situation. The 
following part gives an illustrative example to show the 


thermal performance for such a cascade heat exchange 
network in central heating systems. 


[2+ JC Ce 


(10) 
]C +C. 


[2+ 


G= 


C,+ 1+ 


Illustrative Example and Results 


Case Study of Three Heat Exchangers Network 


In practical central heating systems, the hot fluid tem- 
peratures (This Tho) are often determined by the heat 
source and the overall heating load (Q) is usually known. 
So the main purpose of the cascade heat exchange net- 
work optimization is to find the optimal heating area al- 
locations ((UA)m, (UA), (UA),) as well as the medium 
fluids’ flow rates (Ci, C2), in order to obtain the highest 
cold fluid temperatures (Ti, Te,o) for given conditions. 


Table 1 Known conditions. 
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Table 1 gives the known parameters for a practical 
heating system. Through the aforementioned inverse 
problem and variation method, the optimization results 
can be obtained for this case. Just as Eq. (8) shows, 
without any further constraint condition (e.g., economic 
concern) for these heat exchangers, the heat exchange 
areas should be distributed equally for the three heat ex- 
changers ((UA),=(UA),=(UA),). Table 2 and Fig. 3(a) 
show the temperature variations for each heat exchanger. 


Table 2 Calculation results without economic concern. 


Item Value 
(UA)m, (UA), (UA): (W/°C) 2400 
C (WFC) 710.5 
Cı (WC) 586.9 
Teo CC) 40.1 


Toi CC) 22.1 


It indicates that the medium fluids’ flow rates de- 
creases orderly from the hot to the cold (C.<C\<C,<C;). 
It also shows that under such situation, the inlet and out- 
let temperature differences of the three exchangers are 
equal to each other: 
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Nevertheless, If take the economic factor into consid- 
eration, for instance, suppose the medium heat exchanger 
is much more expensive than the others so that its heat 
exchange area should keep a small and certain value 
((UA)m<(UA), (UA).), the calculation results are shown 
in Table 3 and Fig. 3(b). 


(a) without economic concern (b) with economic concern 
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Fig. 3 Temperature variations in the cascade heat exchange 
network with optimal system parameters. 


Table 3 Calculation results with economic concern. 


Ttem Value 

(UA) m (W/°C) 400 
(UA, (UA): (WPC) 1000 
C (WC) 764.2 

Cı (WC) 554.8 

Teo CC) 33.4 

Toi CC) 15.4 


As Eq. (10) shows, even though the medium heat ex- 
change area should keep a certain value ((UA)m), the oth- 
ers’ heat exchange areas should be also allocated equally 
((UA),;=(UA).>(UA)m). However, compared to the pre- 
vious case results, the cold fluid outlet temperature (7...) 
decreases from 40.1 °C to 33.4 °C and its inlet tempera- 
ture (T,;) decreases from 22.1 °C to 15.4 °C for the same 
known conditions (Q=9000 W, 7,,=70 °C, T,,=60 °C). 
Hence, with the heat exchange area restriction, the ther- 
mal performance of the whole system becomes worse, 
compared to that without such limitation. Moreover, it 
also demonstrates that the mass flow rate difference be- 
tween the two medium fluids increases (C\’<C), C.’>C), 
which makes the temperature differences of the middle 
heat exchanger increase accordingly: 
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Furthermore, Fig. 4 shows that with increasing heat 
exchange area of the middle heat exchanger, the cold 
fluid temperature increases sharply and then decreases 
slightly after the peak value (T, ;=22.1 °C), which can be 
obtained only if the three heat exchange areas are equal 
to each other ((UA),,=(UA),=(UA),). It also indicates that 
the flow rate difference between the two medium fluids 
increases with decreasing heat exchange area of the 
middle heat exchanger. 
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Fig. 4 Parameters variation with changing heat exchange area 
of middle heat exchanger. 


Extension to Multiple Heat Exchangers Network 


The previous case study shows the cascade heat ex- 
change network with three heat exchangers. The prelim- 
inary results indicate that in the ideal situation, the heat 
exchange areas should be distributed equally to each heat 
exchanger and the medium fluids’ flow rate should de- 
crease (Cy>C,) or increase (C,<C,) orderly from the hot 
fluid to the cold one, in order to obtain the highest cold 
fluid temperatures (Toi Teo) for given overall heating 
capacity (Q) and hot fluid temperatures (This Tho). Com- 
paring Eqs. (6) and (9), it can be found that both the two 
cases follow the same principle: 
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Fig. 5 Cascade heat exchange network with multiple heat 
exchangers and medium fluids. 
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Just as Fig. 5 shows, extended to more general situa- 
tion that there are several heat exchangers in series in a 
cascade heat network, in order to obtain the optimal 
thermal performance (i.e., highest cold fluid temperature 
for given hot fluid temperature and overall heating ca- 
pacity), it can be found through the same analysis method 
that the optimal system parameters should meet the fol- 
lowing equation: 


In Eq. (14), (UA); means the heat exchange cost and 


ae represents the fluids flow rate difference 
Ci, C 

of the heat exchanger, respectively. Therefore, to obtain 

the optimal thermal performance of multiple heat ex- 
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changers in series, (UA); E should be equal 
ih 

for each heat exchanger, which also indicates that with 
limited overall heat exchange areas, more heat exchange 
areas (UA); are preferred to be allocated to the heat ex- 
changers where the fluids flow rate difference 
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Conclusions 


In central heating systems, heat is often transferred 
from the heat source to the terminal users by the cascade 
heat exchange network with multi-stage heat exchangers 
in series (heat source, substations and terminal devices). 
In this paper, aim at maximizing the cold fluid tempera- 
tures for given hot fluid temperatures and overall heating 
capacity, the simplified optimization model of the cas- 
cade heat exchange network is established. The optimal 
heat exchange area distribution and the medium fluids’ 
flow rates are determined through inverse problem and 
variation method. The results show that in the ideal situa- 
tion without any extra restrictions, the heat exchange 
areas should be distributed equally for each heat ex- 
changer and the medium fluids’ flow rate should decrease 
(or increase) orderly from the hot fluid to the cold one. It 
also indicates that in order to improve the thermal per- 
formance of the whole system in practical applications, 
more heat exchange areas should be allocated to the heat 
exchanger where the flow rate difference between the 
two fluids is relatively small. 

For the present work, only the simplest case with heat 
exchangers in series is analyzed and discussed to show 
the preliminary application of the inverse problem and 
variation method. In fact, practical heating system is a 
comprehensive heat network with different heat ex- 
changers both in series and in parallel. Additionally, in 
practical engineering fields, not only the fluid flow rate 
itself is of high concern, but also the corresponding pump 
energy consumptions should be taken into consideration 
since the heating users are sometimes far away from the 
heat source. Although the specific results may not be 
applicable to all situations, especially after considering 
more influence factors in real applications, the optimiza- 
tion and analysis approach used here is general, and the 
limitations also raise further study on these topics in the 
future. This work is important for guiding the optimiza- 
tion design of cascade heat exchange networks in central 
heating system. 
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